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Abstract 


The Higgs sector of the minimal supersymmetric standard model (MSSM) 
consists of five physical Higgs bosons, which offer a variety of channels for their 
experimental search. The present study aims to further our understanding of 
the Tevatron reach for MSSM Higgs bosons, addressing relevant theoretical 
issues related to the SUSY parameter space, with special emphasis on the 
radiative corrections to the down-quark and lepton couplings to the Higgs 
bosons for large tan (3. We performed a computation of the signal and back¬ 
grounds for the production processes Wcf) and bbcj) at the upgraded Tevatron, 
with cj) being the neutral MSSM Higgs bosons. Detailed experimental infor¬ 
mation and further higher order calculations are demanded to confirm/refine 
these predictions. 
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I. INTRODUCTION 


The precision electroweak measurements performed at LEP, SLD and the Tevatron are 
consistent with the predictions of the standard model containing a light Higgs boson, with 
mass of the order of the Z boson mass. The searches for such a Higgs particle continue at 
the LEP and the Tevatron colliders. The searches at LEP2 < 200 GeV) are constrained 
by the collider energy, and a Higgs boson with standard model-like properties can be found 
only if its mass is below 105 GeV |l|. 

The potential for discovering a light Higgs boson at the Tevatron collider when it is 
produced in association with a hP or Z gauge boson has been discussed in several studies 
Although the kinematic reach of the Tevatron collider is much greater than for LEP2, 
the backgrounds to Higgs boson searches at hadron colliders are much larger than in e+e“ 
machines. For this reason, large integrated luminosity is essential to establish a signal at 
the Tevatron. Within the standard model, the general conclusion is that Run H, with a 
total integrated luminosity of about 2 fb“^ per detector, will be unable to extend the Higgs 
boson mass reach of LEP2. The main questions are: what is the theoretical motivation for a 
Higgs boson with a mass slightly above the LEP2 reach, and what is the necessary upgrade 
in luminosity to cover that region? 

We address the theoretical motivation by appealing to the minimal supersymmetric ex¬ 
tension of the standard model (MSSM). The MSSM has the remarkable property that, for a 
sufficiently heavy supersymmetric spectrum, it fits to the precision electroweak observables 
as well as the standard model 0. Moreover, the lightest GP-even Higgs boson mass mh is 
constrained to satisfy ra^, ^ 130 GeV |§,^. The Higgs sector of this model consists of two 
Higgs doublets, with two GP-even Higgs bosons, h and if, one GP-odd Higgs boson. A, 
and one charged Higgs boson, H^. This richer spectrum allows for different production and 
decay processes at LEP and the Tevatron colliders than in the standard model. 

In the supersymmetric limit, the neutral components of the two Higgs boson doublets 
Hi and H 2 couple to down- and up-type quarks, respectively. Lepton fields couple only to 
the Hi Higgs boson. The MSSM, tree-level Yukawa couplings of the down quarks, leptons 
and up quarks are related to their respective running masses by 


hd ~ 


rrid 

V cos f3 ’ 


hi ~ 


mi 


V cos [3 ’ 


, mu 

h ~ - 

i;sin p 


( 1 . 1 ) 


where tan/? = V 2 /V 1 is the ratio of the vacuum expectation values of the two Higgs doublets, 
and V = \Jvi-\- ^1=174 GeV. In the standard model, only the top quark Yukawa coupling 
ht is of order one at the weak scale. In the MSSM, instead, the bottom and r Yukawa 
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couplings, hh and hr, can become of the same order as ht, if tan/? is sufficiently large. This 
can have important phenomenological consequences. 

Quite generally, the two CP-even Higgs boson eigenstates are a mixture of the real, 
neutral Hi and H 2 components, 

h \ I — sin a cos a \ / H^ 

H I \ cos a sina } \ iP® 



and the lightest CP-even Higgs boson couples to down quarks (leptons), and up quarks by its 
standard model valnes times — sina/ cos/5 and cos a/sin/?, respectively. The conplings to 
the heavier CP-even Higgs boson are given by the standard model valnes times cos a/ cos (3 
and sin a/ sin (3, respectively. Analogously, the coupling of the CP-odd Higgs boson to down 
quarks (leptons) and np qnarks is given by the standard model conpling times tan f3 and 
1/tan/?, respectively. Moreover, the lightest (heaviest) CP-even Higgs boson has ZZh and 
WWh [ZZH and WWH) conplings which are given by the Standard Model value times 
sin(/5 — a) (cos(/5 — a)), while it can be prodnced in association with a CP-odd Higgs boson 
with a ZhA (ZHA) conpling which is proportional to cos(/5 — a) (sin(/5 — a)). 

For snfficiently large valnes of the CP-odd Higgs boson mass tha, the effective theory at 
low energies contains only one Higgs donblet, with standard model-like properties, in the 
combination 


h ~ = FT/ cos (3 + Hi sin (3- < 0^^ >= v , (1.3) 

where sina ~ —cos {3, cos a ~ sin/? and hence sin^(/5 — a) = 1. In this limit, if all 
supersymmetric particles are heavy, all phenomenological conclusions drawn for a SM Higgs 
boson are robust when extended to the the lightest CP-even Higgs boson of the MSSM, 
which, as mentioned before, is at the same time constrained to have a mass below abont 
130 GeV. For large tan (3, one of the CP-even neutral Higgs bosons tends to be degenerate 
in mass with the CP-odd Higgs boson and conples strongly to the bottom qnark and tan 
lepton. The other CP-even Higgs boson has standard model-like conplings to the gange 
bosons, while its conpling to the down qnarks and leptons may be highly non-standard. If 
mA is large, then h is the Higgs boson with SM-like properties as described above. If ruA is 
small, then H is the one with SM-like conplings to the gange bosons. In the following, the 
symbol 0 denotes a generic Higgs boson. 

In this article, we analyze the discovery potential of the Tevatron collider for MSSM 
Higgs bosons in different production channels. Section 2 discusses the signals from Higgs 
boson prodnction in association with weak gange bosons and their dependence on the MSSM 
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parameter space. Section 3 contains details of Yukawa coupling effects for large tan /3. 
Section 4 deals with the phenomenological implications of these effects for signals in the W(f) 
and bb(j) production channels. In section 5, we consider the correlation between the bottom 
mass corrections and the supersymmetric contributions to the branching ratio B{b —>■ sy). 
Finally, section 6 is reserved for our conclusions. 

II. SIGNALS FROM W<p AND Zcp PRODUCTION 

The production of hF0, followed by the decays W{—> ez/g) or IF(—>• pz/^) and bb), is 
the gold-plated search mode for the standard model Higgs boson (j) = cf)^^ at the Tevatron 
collider, while LEP2 is sensitive to the process. The production of Z(j) may also 

be useful at the Tevatron, depending on the efficiency for triggering on missing Et {^t) 
and the bh mass resolution. Additionally, the all-hadronic decays of hF0 + Z(p may extend 
the reach, and the Higgs boson decay (p{^ r+r“) may be observable. However, there are 
several, unresolved experimental issues concerning these channels that require detailed study 
by the experimental collaborations. For this reason, at present we shall only consider the 
Wcj) channel at the Tevatron. 

To quantify the experimental reach in a model independent way, it is useful to consider 
the function 

ajpp^W^) B{(t)^bb)^ 

a{pp B{(i)SM ^ 55) ’ ^ ^ 

where a denotes a production cross section, B denotes a branching ratio and (j)^^ represents 
a standard model Higgs boson. In the MSSM, the ratio of cross sections is just given by 
sin^(/3 — a) or cos^(/3 — a) depending on 0 being the lightest or heaviest CP-even Higgs, 
respectively, Q while the ratio of branching ratios has a more complicated behavior. It 
is important to notice that, barring the possibility of large next-to-leading-order (NLO), 
SUSY corrections to the WW(f) vertex, there is no enhancement of the production cross 
section in the MSSM over the standard model. On the other hand, the branching ratio 
to bb and final states are affected by the factors — sin a/ cos (3 for h and cos a/ cos {3 

for H over the SM Higgs boson couplings to down quarks and leptons. These factors can 
produce an increase or decrease of the MSSM coupling of the Higgs boson to bottom quarks. 


^Observe that sin(/3 — a) (cos(/3 — a)) denotes the component of the lightest (heaviest) CP-even 
Higgs boson in the combination which acquires a vacuum expectation value, Eq. O- 
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depending on the value of the CP-odd mass, tan (5 and the top and bottom squark mass 
parameters. In this study, the Higgs boson properties are calculated using the program 
Hdecay |TO . 


As mentioned above, for large m^, the low-energy, effective theory contains only one 
Higgs boson with SM-like properties. The Yukawa couplings tend to the standard model 
values, and sin^(/? — a) ~ 1. As long as no new decay modes are open, h has the same 
properties as 0'^^, and a discovery or exclusion limit for a standard model Higgs boson 
applies equally well to h, and R = 1. If Higgs boson decays to sparticles become important, 
then it is quite likely that additional Higgs boson production modes exist and enhance the 
potential signal, rather than decrease it. For example, if the sparticle spectrum is of the order 
of nih, then processes like pp —> A" 2 (—> Nih)Ci{^ NiW) + X can occur. In our analysis, 
we shall always consider the limit of heavy sparticle masses, where such supersymmetric 
contributions to the Higgs production and decay processes are negligible. 

In the large ttia limit, the renormalization-group improved result for the lightest Higgs 
boson mass, including two-loop leading-log effects |]8|Jp, pT| , p!^ , has the approximate analytic 
form [R|: 


ml = mJsM 


= Ml coi? {213) (1 - 
3 


4 

rui 


4^2 ^2 
Xf = 2a^ 


1 


t 






(Xtt + 


1 - 


a 

12 , 


a = At — p/ tan/?. 


( 2 . 2 ) 


where p = p/Ms and At = At/Ms, with M| = (m|^ + mU/2, and t = \og{Mg/ml). The 

above formula is based on an approximation in which the right-handed and left-handed 

stop supersymmetry breaking parameters are assumed to be close to each other, and hence 

the stop mass splitting is induced by the mixing parameter m^a x Ms- Moreover, this 

|m? —I 

expression is based on an expansion in powers of mta/Ms and is valid only if ^2 , < 0.5, 

where and are the lightest and heaviest stop mass eigenstates. This simplihed 
expression is very useful for understanding the results of this work, although we go beyond 
this approximation in our full analysis. Finally, in the above, we have ignored corrections 
induced by the sbottom sector, which, as we shall discuss below, may become relevant for 
very large values of tan /?. 

The value of mh in Eq. (2.2) is maximized for large values of tan/3 and Ms and = 6. 
Due to the dependence of the lightest CP-even Higgs mass on tan /3, LEP will be able 
to probe the low tan (3 region of the MSSM. Indeed, recent analyses suggest that even the 
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present, relatively low bounds on a standard model-like Higgs boson from LEP2 have strong 
implications for the minimal supergravity model |^. Moreover, it has been shown that, in 
the large region, LEP2 will probe values of tan/3 < 2 for arbitrary values of the stop 
masses and mixing angles [Q. Since, in general, the lower bound on tan/3 is expected to be 
obtained for large values of ttla, the absence of a Higgs signal in the ZZh channel at LEP2 
will provide a strong motivation for models with moderate or large values of tan /3. 

For lower values of the CP-odd Higgs mass, sin^(/3 — a) can take any value between 0 
and 1, and it is a model-dependent question whether h or H produces a viable signal in the 
W(j) production channel. For moderate or large values of tan/3 5, it is easy to identify the 
main properties of the CP-even Higgs sector. More specihcally, three cases may occur: 


a) If niA < m^sM, then sina ~ —1, cosa ~ 0[1/ tan/3) and cos(/3 — a) ~ 1. In this case, 
the heaviest CP-even Higgs boson has a production rate which is similar to the standard 
model case. The branching ratio of the decay into bottom quarks and r leptons, however, 
can become highly non-standard, since cosa and cos/3, may differ by a factor of order one. 


b) If niA > m^sM, then the lightest CP-even Higgs boson has a production rate similar to 
the standard model case. For uia ^ mfpsM, the branching ratio of the decay of this Higgs 
boson into down quarks is standard model-like. However, when becomes close to m^sM, 
there can be important differences in the branching ratios with respect to the SM ones. 


c) If rriA — m^sM, then sin^(/3 — a) ~ cos^(/3 — a) ~ 0(0.5), and the couplings of both 
neutral CP-even Higgs bosons to bottom quarks tend to be highly non-standard. 

To better understand the behavior of the Higgs boson branching ratios in these different 
cases, we analyze the Higgs boson mass matrix. Assuming the approximate conservation 
of CP in the Higgs sector, the CP-even Higgs masses may be determined by diagonalizing 
the 2x2 symmetric mass matrix After including the dominant one-loop corrections 
induced by the stop and sbottom sectors, together with the two-loop, leading-logarithm 
effects, the elements of are [§ 


~ m\ /3 + M| cos^ /3 


’H ^- 2*2 n~2 

Sin Pa 


1 + 


IbTT^ 


(eh* - 2hl - IQg. 


+ 0{hlMl 


AI 22 — cos^ /3 + M| sin^ /3 

(l.bh] + Q.bhl - 8gl) 




1 + 


Ibvr^ 
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+ Attt ( 1 — 


12 , 


1 + 


IGvr^ 


(3/i? + hl- 16gl 


Aht 


IGtt^ 


sin^ Pfj. 


1 + 


(9hl-5h^,-16gl) 




m\ + M| 


hy 

Stt^ 


IGtt^ 

{y - 


+ oyMi) 


sin f3 cos f3 


h‘^v‘^ „ r - ~ 

sin^ f3jia Ata - 6 

lOTT"' L 


+ 


yui 

327r2 


/iO 


1 + 


167r2 


(4.5/1? - 0.5/1? - 16g 


(2.3) 


where g^ is the QCD running coupling constant. In the above, we have assumed, for sim¬ 
plicity, that Afc ~ 0 (an assumption we shall always make in the following analysis) and 
retained only the leading terms in powers of and tan f3. We have also included the 
small, OyM"^) correction to explicitly because it plays a relevant role in our anal¬ 
ysis. The above expressions hold only in the limit of small splittings between the running 
stop masses. Moreover, the condition 2mt ma.x{\At\,\fi\) < M| must be fulhlled. Similar 
conditions should be fulhlled in the sbottom sector. The leading, two-loop, logarithmic 
corrections to the squared Higgs mass matrix elements included above can be as large as 
10 — 20 % when supersymmetric particles are heavy, and are very relevant in determining 
the Higgs boson mass eigenvalues and mixing angles. Observe that Eq. (|2.2| ) may be easily 
obtained from the above expression, by computing the determinant of the Higgs boson mass 
matrix and setting the heavy CP-even Higgs boson mass approximately equal to mA^ 

The mixing angle a can be determined from the expression 

sin a cos Q = —^ =. (2.4) 

In the limit that AI 12 —>■ 0, either sin a or cos a ^ 0. For moderate or large values of tan/?, 
if case a) is realized and cos a ~ 0, the coupling of the standard model-like Higgs boson| 
H to bb and is diminished, and decays to gg, 77 , cc, and can be greatly 

enhanced over standard model expectations [^. The same can happen for h in case b) 
when sin a ~ 0. For moderate or large values of tan/?, the vanishing of M .\2 leads to the 
approximate numerical relation: 


^There is, however, a slight discrepancy in the subdominant, Yukawa-dependent, two-loop, 
leading-logarithmic corrections, which is due to the fact that the expressions written above are 
only strictly valid for values of the CP-odd mass of the order of the weak scale. 

^From now on, the term “standard model-like Higgs boson” refers to a Higgs boson with standard 
model couplings to the gauge bosons with no implication about its couplings to fermions. 
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ml 


27r^ 


(s/i^ - 


+1 


tan/? |- 
100 


fia{2aAt — 11 ) 


1 - 


15 

IOtt^ 


(2.5) 


where we have neglected the bottom Ynkawa conpling effects and replaced and and 
the weak gange conplings by their approximate nnmerical valnes at the weak scale. For low 
valnes of or large valnes of the mixing parameters, a cancellation can easily take place 
for large valnes of tan/3. For instance, if Ms — 1 TeV, /2 = —d = 1, and — 80 GeV, a 
cancellation can take place for tan (3 ~ 28, with the spectrnm ~ and mu — 117 GeV. 
The heaviest GP-even Higgs boson has standard model-like conplings to the gange bosons 
(cos^(/? — a) ~ 1 ), bnt the branching ratios for decays into bosons, glnons and charm 
qnarks are enhanced with respect to the SM case: Bw = -34, = .27 and Be = .11. For 

the same valne of bnt larger valnes of the stop mixing parameters, ft = d = a/7 (at the 
edge of the region of validity of the above approximation), an approximate cancellation of 
the tree-level bottom and r lepton conplings is achieved for tan /3 = 20, for which mn — 124 
GeV, with branching ratios B]y = .57, Bg = . 21 , and Be = .08. 

An interesting point is that, for large valnes of tan/? and valnes of the stop mixing 
parameters which maximize A422, (d = VG), the dominant, m^-dependent corrections to 
the off-diagonal elements of the Higgs boson matrix vanish. In this case, the corrections to 
dominated by the dependent terms (see Eq. (^)), which, for Ms — 1 TeV, 

cannot be large enongh to indnee an approximate cancellation of the off-diagonal terms for 
\fi\ < 1 in the region of mA ~ 400 GeV and tan/? < 50 considered in this article. However, 
even for \fl\ = 1, the impact of the radiative corrections to the off-diagonal elements of the 
Higgs boson mass matrix may be very relevant for low valnes of mA, and we expect large 
variations of the branching ratio of the decay of the heaviest GP-even Higgs boson into 
bottom qnarks with respect to the choice of the sign of /i in this region of parameters. 

Moreover, away from |a| = a/ 6 , the m^-dependent radiative corrections to AI 12 depend 
strongly on the sign of fl x At [At ~ d for large tan /? and moderate /i) and on the valne of 
\At\. For the same valne of a, a change in the sign of /i can lead to observable variations 
in the branching ratio for the Higgs boson decay into bottom qnarks. If \At\ ^ y^ 11/2, the 
absolnte valne of the off-diagonal matrix element, and hence, the conpling of bottom qnarks 
to the standard model-like Higgs boson tends to be snppressed (enhanced) for valnes of 
/i X < 0 (/t X > 0). For larger valnes of \At\, instead, the snppression (enhancement) 
ocenrs for the opposite sign of /Z x A^. Finally, it is important to stress that, in the large 
tan /? regime, extra corrections to the Ynkawa conplings may be important depending on 
the MSSM spectrnm, and we shall come back to this topic later in Section 3. 
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A. W(^ bb) signal and backgrounds 


The signal W (—> >■ 66 ) contains two real 6 -jets, which can be used to distinguish 

it from many potential backgrounds. A remarkable increase in the double tag efficiency over 
the Run I estimate is expected to be accomplished by loosening the requirements for the 
second tag [§]. This is possible since the first 6 -tag already significantly reduces the fake 6 
background. In the numerical analysis we assume a double tagging efficiency for Run II of 
= .52 X 1.8 = .45. 

Given the expected high 6 -tagging efficiency and the low mistagging rate, the most 
important backgrounds are those with a real charged lepton and two real 6 taggable jets. 
The backgrounds considered are hh),W^Z^{-^ bb),t{—>- bW~^)t{—>- bW~),W~^* —> 

6 IT+) 6 , and qg —> qt{-^ 6 IR’'') 6 , and hermitian conjugate (h.c.) processes where 
appropriate. The boson decays leptonically, where or /i^, except 

in the ti background, where the second boson can decay hadronically or to a r. 

Events are required to have one lepton with p^rp > 20 GeV and < 2. The lepton 
must be isolated from jets with a separation = v^A0^”+~A^ > .4, where A0 and Ag 
are the difference in azimuthal angle and pseudorapidity between the lepton and jets. Jets 
which can be heavy flavor tagged must have > 15 GeV and < 2. Additional jets 
are resolved if > 15 GeV and < 2.5, and leptons if pP > 10 GeV and < 2. 
Jets must be separated from each other by > .7. A Higgs boson signal is defined as an 
excess above backgrounds in the invariant mass distribution of the bb pair, A mass 

resolution of about aM = is assumed. For Mh = 100 GeV, this means a window 


84 < < 116 GeV. 

To further enhance the signal over the background, we apply additional cuts. The angle 
between the and h in the W^h center of mass system, 9h, can be exploited |T^. The 
dominant background from bb) tends to peak at cos 6h = ± 1 , and a cut of | cos 9h\ < 

.8 is optimal. Top quark pair production events, which are a large potential background, 
produce an additional boson, which can decay hadronically or leptonically. The extra 
decay products from this W decay can be used to veto such events. Vetoing events with at 
least one jet with Et > 30 GeV, a pair of jets separately having Et > 15 GeV, or extra 
leptons with px > 10 GeV successfully limits this background. 

The estimate of the signal and background used in this analysis are based on a parton 


level calculation |]^, and the final state partons are interfaced to a detector simulation 
to account for finite detector resolution in measuring energies and angles |]T^. The NLO 
QGD corrections to those processes which are order at tree level are large. The order 
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correction to production at Tevatron energies is about 1.3 |]T^. A similar enhancement 

occurs for the signal process which has the same initial state, and the rate for the 

process ud —> > tb, with both initial and final state corrections, is enhanced by 

a factor of 1.7 |2^ over the lowest result using CTEQ3L [^] structure functions. The 
single top production processes qb q't increases by a factor of about 1.3 Results are 
normalized to these numbers. In addition, ti and W^g*{^ bh) are evaluated at the scale m^, 
which gives good agreement with the present data. However, higher-order calculations of 
the W^g*{^ bb) production rate and kinematics are necessary to conhrm our understanding 
of the standard model backgrounds p3 . 


The results of this analysis will be used below. These results are in good agreement with 
other studies after accounting for different assumptions for the mass resolutions au 
(see also Ref. 


B. Results on Wh and WH 

The signal and background estimates from the analysis described above can be used to 
estimate the discovery or exclusion reach of the Tevatron. For a hxed integrated luminosity 
and a Higgs boson mass, we can determine which values of R would lead to a discovery 
with a 5cr signihcance (a = S/\fB) or a 95% C.L. (1.96 ct) exclusion. If the R-contour lies 
below R = 1, then a standard model-like Higgs boson could be discovered or excluded. The 
exclusion potential of the W(f) channel at the Tevatron is summarized in Fig. |^, while the 
analogous discovery potential is described in Fig. |^. Fig. shows the 95% C.L. exclusion 
limit as a function of for LFP2 running at ^/s = 192 GeV and collecting 150 pb“^ of 
data (dash-dot) and for the Tevatron with 30 fb“^ (solid), 10 fb“^ (long-dash) and 2 fb“^ 
(short-dash). 2 The sensitivity of the numerical results for the Tevatron are demonstrated 
by the markers, which show the 95% C.L. exclusion limit for 2 fb“^ from a different study 
that uses the Run I CDF mass resolution (squares) and from the same study with the 
backgrounds doubled (triangles). Fig. || shows 5a discovery curves as a function of 
for LFP2 and the Tevatron under the same assumptions as above. Assuming improved 
mass resolution from Run I, the Tevatron Run H, with 2 fb“^ of data, may exclude a 110 


^This gives a conservative estimate of the LEP2 reach, which is expected to be better: for ^/s = 200 
GeV and 200 pb“^, a maximal discovery reach of 105 GeV is expected for a standard model-like 
Higgs boson. 
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GeV standard model Higgs boson at the 95% C.L. However, assnming the present mass 
resolntion, the exclnsion limits for 2 fb“^ can range from abont 90 to 102 GeV, assnming 
the backgronnd and experimental efficiencies are nnderstood within a factor of 2 in the 
present stndies. With 30 fb“^ and improved mass resolntion, a SM Higgs boson of mass 130 
GeV can be discovered within this channel. In order to analyze the MSSM case, we shall 
consider the case of improved mass resolntion, since, as we shall show, this will be necessary 
to have good coverage of the MSSM Higgs sector. 

The i?-contonr as a fnnction of Higgs boson mass is the starting point also for a MSSM 
Higgs boson analysis. After specifying the parameters of the top sqnarks, the fnnction R 
can be calcnlated as we scan throngh uia and tan (3. For the MSSM, the nnmerical resnlts 
are illnstrated in Figs. Fignres correspond to three common choices of the MSSM 
parameters. This set is not exhanstive, bnt is meant to illnstrate the effect of the stop 
mixing parameters when all sparticles are relatively heavy. We have taken the sqnarks to 
have masses Ms = 1 TeV. The higgsino mass parameter is taken to have the valnes ±1 TeV. 
Finally, the stop trilinear conpling At is chosen so that the stop mixing parameter is either 
very small (minimal mixing), d ~ 0, or, in the limit of large m^, it maximize the lightest 
GP-even Higgs mass (maximal mixing), |d| ~ ^/6. 

In the case of maximal mixing, the radiative corrections to the Ml 12 depend on the sign 
oi ft X d. This dependence is obvions in Figs. and where the discovery reach for the 
cases /r = ±1 TeV and d = a/6 is displayed. An increase or decrease of the effective Ynkawa 
conpling of a standard model-like Higgs boson to bottom qnarks, for negative or positive 
valnes of /r x d, leads to large variations in the branching ratio to b qnarks. This has an 
important impact on the Inminosity reqnired to observe a Higgs boson in this channel. For 
Higgs boson masses in the range 125-130 GeV, as expected for the SM-like Higgs boson in 
the limit of large tan (3 and maximal mixing, a decrease in the branching ratio to bottom 
qnarks is compensated by an increase in the branching ratio to W^W^, 

From Figs. ^ and ^ it follows that the dependence on the sign of /i x d is particnlarly 
strong in the regime of low valnes of m^, bnt, as we discnssed before, might become relevant 
even for relatively large valnes of ttia — 0(300 GeV). The snppression of the conpling which 
is obtained for low valnes of tua and large valnes of tan f3 leads to a problem for detecting 
the heavy Higgs boson in this regime for positive valnes of /i x d (see Eq. (0)). 

The region of — 120 GeV {niA — m^sM) and large tan/? is difficult to observe, since 
this is the region of maximal mixing and sin^(/? — a) ~ cos^(/? — a) ~ 0(0.5). Although 
this limitation can be overcome with larger luminosity, a window of non-observability would 
remain for both signs of /i x d. Fortunately, if tan (3 is sufficiently large, the two GP-even 


11 


Higgs bosons tend to have similar masses. Whenever the mass difference is less than 10 
GeV, we have increased the mass window to include both signals. Usually, the number of 
events within ±2 (Tm of a given Higgs boson mass is used to quantify the significance of 
any deviation from the expected mean number of events. To combine the signals from two 
Higgs bosons with slightly different masses, the mass window extends from — 2aM to 
+ 2aM, where (pio and (phi are the lighter and heavier of the two Higgs bosons. Using 
this procedure, we can extend the coverage for very large tan/? and tua — m^sM. However, 
a window of non-observability remains for tan/? ~ 5, since the mass difference is large and 
sin(/? — a) is suppressed in this region. | Finally, the region of ijia — 150 GeV is easier to 
observe, since sin^(/3 — a) is already of order one in this region and the bottom coupling of 
the lightest GP-even Higgs boson is strongly enhanced with respect to the standard model 
case, implying an increase in the branching ratio of this Higgs boson into bottom quarks. 

In Fig. the case of minimal mixing (a = 0, /2 = 1) is displayed. In the large ttia 
limit, the lightest GP-even Higgs mass is of order 110-115 GeV for moderate or large values 
of tan/? and hence detectable for luminosities of order 10 fb“^. The characteristics of this 
case are similar to the case of maximal mixing, although, due to the lower values of the 
lightest GP-even Higgs mass, lower luminosities are required to cover the large niA region 
and, in addition, the window of non-observability for 32 fb“^ shrinks to a very small region 
of parameter space for ~ m^sM. For minimal mixing, the results are insensitive to the 
sign of /i. This occurs since the dependence of the Higgs boson mass matrix (and hence of 
the GP-even Higgs masses and their couplings to bottom quarks) on the sign of fi arises 
through the radiative corrections to the off-diagonal elements, which are proportional to 
fl X d and vanish when a ~ 0. 

To illustrate the sensitivity of our results to experimental resolution, we have constructed 
Fig. H from the results of another study ^ which used the present GDF mass resolution. For 
worse mass resolution, the discovery reach at large niA and near ttia — m^sM is compromised, 
but the general features are the same. 

In Figs. 1^ and ||, we present cases in which the stop mixing parameters are such that 
the bottom Yukawa coupling of the standard model-like Higgs boson can be efficiently 
suppressed in a large region of parameters. For this, we have taken values of the mixing 
parameters At = —/x = 1 and 1.5 TeV. In these cases, in the limit of large values of tua and 


®If we had not combined the signals, the region of non-observability for tha — rn^sM would have 
extended to tan/? ~ 50 (30) in Fig. 3 (4). 
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for moderate or large values of tan/3, the lightest CP-even Higgs boson mass takes values 
in the range 115-120 GeV and 120-125 GeV, respectively. In these two cases, windows of 
non-observability appear associated with the suppression of the bottom Yukawa coupling of 
the standard model-like Higgs boson, i.e. vanishing A4i2- For /i = 1 TeV, the mass of the 
standard model-like Higgs boson is smaller and the bottom Yukawa coupling cancellation 
is more difficult than in the case of /r = 1.5 TeV (see Eq. (|2.5|) ). Hence, although in the 
former case the window of non-observability is small and restricted to small values of uia, 
in the latter case, for large values of tan/3 the window of non-observability extends up to 
relatively large values of the CP-odd Higgs mass. It would be very interesting to check 
if these windows may be efficiently covered by using the IVIV decay mode of the Higgs 
boson |T3| , ^ . 

It is clear from Figs. ^|-§that a single detector at the Tevatron requires about 30 fb“^ for 


a reasonable coverage of the MSSM parameter space, far beyond the region already covered 
by LEP2. Other decay channels besides h,H bb need to be explored to cover specihc 
regions of parameter space. For values of Ms — 1 TeV and moderate or large values of 
tan (3, the standard model-like Higgs boson tends to be heavier than 100 GeV. In this case, 
the LEP reach in the Zh and ZH channels is highly reduced, and most of the coverage 
usually shown is induced by the hA production (since HA is kinematically limited) .0 An 
essential advantage of the Tevatron is the fact that it can overcome this kinematic limitation 
and give a signihcant coverage of the m^-tan/3 plane via the Wh and WH channels even 
for large values of tan/3. Obviously, the addition of the Zh and ZH channels would be 
useful to conhrm a signal, or to enhance the possibility of a discovery with lower integrated 
luminosity. 


III. YUKAWA COUPLING EFFECTS IN THE LARGE tan /3 REGIME. 

In the SM, the coupling of the Higgs boson to b quarks is proportional to the bottom 
Yukawa coupling = rrib/v. Within the MSSM, the effective bottom Yukawa coupling can 
be quite different than in the standard model case. This is due not only to the dependence 
on the Higgs mixing angles, discussed above, but also to the presence of large radiative 
corrections in the coupling of bottom quarks and r leptons to the neutral components of the 


®For yG = 200 GeV and 200 pb LEP2 can discover a Higgs boson in the hA channel for 
vriA < 90 GeV and large tan j3. 
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Higgs doublets, that lead to modifications of the relation, Eq. ( p..l|) , between the bottom 
Yukawa coupling and the running bottom mass p6|-p8|. To better understand this, it is 
necessary to concentrate on the properties of the large tan (3 regime. In this regime, to 
a first approximation, only H 2 acquires a vacuum expectation value (here H 2 denotes the 
whole SU{2)l doublet). This means that this Higgs doublet contains the three Goldstone 
bosons and a neutral Higgs boson with standard model like couplings to the electroweak 
gauge bosons. The other Higgs doublet Hi does not communicate with the electroweak 
symmetry breaking sector and contains, in a first approximation, a CP-even and a CP-odd 
Higgs field, which are almost degenerate in mass, and a charged Higgs field, whose mass 
differs from niA only in a D-term, = m\ + M^. The SM-like Higgs boson acquires a 
mass given by 


Tn^sM — MI 2 ~ M| + radiative corrections. 


(3.1) 


and its dependence on is suppressed by a 1/ tan^/3 factor (see Eqs. (p.2|) and ( p.3|) ). 

Supersymmetric one-loop corrections to the tree-level, running bottom quark mass can 
be significant for large values of tan (3 and translate directly into a redefinition of the relation 
between the bottom Yukawa coupling entering in the production and decay processes and the 
physical (pole) bottom mass. Some of the phenomenological implications of these corrections 
have been considered for MSSM Higgs boson decays [^. The main reason why these one- 
loop corrections are particularly important is that they do not decouple in the limit of a 
heavy supersymmetric spectrum. As mentioned above, in the supersymmetric limit, bottom 
quarks only couple to the neutral Higgs Hi. However, supersymmetry is broken and the 
bottom quark will receive a small coupling to the Higgs H 2 from radiative corrections. 


C ~ hH^bb + AhH^bb. 


(3.2) 


The coupling Aht, is suppressed by a small loop factor compared to hb and hence, one would 
be inclined to neglect it.[] However, once the Higgs doublet acquires a vacuum expectation 
value, the running bottom mass receives contributions proportional to AhbV 2 . Although Ahb 
is one-loop suppressed with respect to hb, for sufficiently large values of tan [3 {y 2 3> Ui) the 
contribution to the bottom quark mass of both terms in Eq. (|3.2|) may be comparable in 
size. This induces a large modification in the tree level relation, Eq. (0), 

mb = hbVi{l +A{mb)), (3.3) 


^ In the above we are explicitly neglecting corrections to the hb coupling of 0{Ahb). 
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where A(mfc) = Ahbtan/3/hfc. 

The function A{mf,) contains two main contributions, one from a bottom squark-gluino 
loop (depending on the two bottom squark masses and and the gluino mass Mg) 
and another one from a top squark-higgsino loop (depending on the two top squark masses 
and and the higgsino mass parameter fi). The explicit form of A{mb) at one- 
loop can be approximated by computing the supersymmetric loop diagrams at zero external 
momentum [Ms 3> mb) and is given by 


2q/o y, 

A{mb) ~ —Mg/itan/? /(M^^, Mg) + —A^/itan/? fi), 


Svr 


(3.4) 


where 0^3 


^|/4vr, Yt 


hi 

An ’ 


and the function / is given by, 


I{a,b, c) 


alb'^ ln(a^/6^) + ln(6^/c^) + ln(c^/a^) 

{a? — b‘^){b‘^ — c^){a? — c^) 


(3.5) 


and is positive by dehnition. Smaller contributions to A{mb) have been neglected for the 
purpose of this discussion. It is important to remark that these effects are just a manifes¬ 
tation of the lack of supersymmetry in the low energy theory and, hence, A{mb) does not 
decouple in the limit of large values of the supersymmetry breaking masses. Indeed, if all 
supersymmetry breaking parameters (and /i) are scaled by a common factor, the correction 
A(mfe) remains constant. 

Similarly to the bottom case, the relation between m,- and the r lepton Yukawa coupling 
hr is modihed: 

nir = hrVi{l + A{mr)). (3.6) 

The function A{mr) contains a contribution from a tau slepton- bino loop (depending on 
the two stau masses Mf^ and Mr^ and the bino mass parameter Mi) and a tau sneutrino- 
chargino loop (depending on the tau sneutrino mass Mj>^, the wino mass parameter M 2 and 
fi). It is given by the expression 

A(mT-) = —Mi/itan/?/(M~ , M™, Ml)—M 2 /itan/?/(My ,M 2 ,p.), (3.7) 

dvr dvr 

2 

where ai = gi is the U{1) hypercharge coupling, a ;2 = g 2 is the SU{2) weak isosopin 
coupling. 

Since corrections to hr are proportional to ai and 0 ) 2 , they are expected to be smaller 
than the corrections to hb- Although the precise values of A{mb) and A{mr) are model 


27,28 
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dependent, the leading term in the tan mass corrections has a factor M 2 ,/i) oc 

M 2 /i/(max(m?^, M|,/i^) < 1, and hence, for tan/3 < 50, /S.{mr) < 0.15. In the following, 
we consider the impact of the bottom mass corrections assuming A{mT-) -C A(mfe), and 
using the expression A{mb) = Ktan(3 to parametrize possible radiative corrections. Since 
the value of at the scale Ms is of order 0.1, and if all soft supersymmetry breaking 
parameters and /i are of order of 1 TeV, the coefficient K can have either sign and will be 
of order \K\ ~ 0.01. One can also consider cases in which the bottom mass corrections are 
highly suppressed. This happens naturally in the case of approximate R and Peccei-Quinn 
symmetries in the theory, which make the value of the gaugino masses and the stop mixing 
parameters much lower than Ms . 


It is instructive to return to the couplings of the lightest and heaviest CP-even Higgs 
bosons and of the CP-odd Higgs boson to bottom quarks. The CP-odd Higgs boson coupling 
to bottom quarks is given by 


c = Ab-nb 


(3,8) 


with 

T77 u 

= hb sin (3 + Ahb cos /3 ~ h;, sin = -——— tan j3. (3.9) 

(1 + A[mb))v 

Using Eqs. (|3.2|) and ( |1.2| ), together with the relation of the bottom Yukawa coupling to 
the bottom mass, Eq. ( p.3|) , it is easy to show that the effective couplings of the CP-even 
Higgs bosons, hb and hb, 

C = hbhhh + hbbbH (3.10) 

are approximately given by 

- mb sin a 

hb ^ - - 

V cos p 


A{mb) 

1 + A{mb) 


1 + 


tan a tan (3 


(3.11) 


~ mb cos a 
hb ^^ 
V cos (3 


A(mfc) 

1 + A(mfe) 



(3.12) 


The value of A{mb) in eq. (|3.4|) is dehned at the scale Ms, where the sparticles are decoupled. 
The hb and Ahb couplings should be computed at that scale, and run down with their 
respective renormalization group equations to the scale m^, where the relations between 
the couplings of the bottom quark to the neutral Higgs bosons and the running bottom 
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quark mass, Eqs. (|^) , (|3.11|) , and (|3.12|) are defined. In the present study we have 
dehned the running bottom mass at the scale ttia as a function of mh{mh) — 4.25 GeV while 
using two-loop renormalization group equations in the effective standard model theory at 
scales Q < rriA- The above procedure leads to a consistent dehnition of the bottom quark 
couplings to the Higgs bosons when the three neutral Higgs boson masses are of the same 
order. For large values of the CP-odd Higgs boson mass, instead, hj, must be evolved with 
SM renormalization group equations from niA down to rrih- The dehnition of the couplings 
of the bottom quark to the Higgs bosons at the scale of the corresponding Higgs boson mass 
is chosen to take into account the bulk of the QCD correction. Indeed, it is known that this 
choice of scale represents well the bulk of the QCD corrections to the Higgs boson decay 
into quarks and gluons [|^. However, for the production process bbcp, a complete study of 
the NLO effects remains necessary to get a dehnitive estimate of the Tevatron reach in this 
production channel. 

It is interesting to study different limits of the above couplings. For large tua 3> m^sM, 
the lightest CP-even Higgs boson should behave like the SM particle. This is fulhlled since, 
in this limit cos/? ~ — sin a and sin/? ~ cos a. Hence, hb = fUb/v, which is the standard 
coupling. In the same limit, hb — hbsin/?(l + 0{A{mb)/ tan^/?)). Even in the presence 
of radiative corrections to the bottom quark couplings, the heaviest CP-even Higgs boson 
coupling is approximately equal to the CP-odd one. When starts approaching m^sM, the 
above relations between the angles a and /? are slightly violated. Due to the large 1/cos/? 
factor appearing in the dehnition of the Yukawa coupling, Eq. ( ^.11| ), a small departure 
from the above relations can induce large departures of the coupling hb with respect to the 
standard model value. For -C m^sM, instead, sin/? ~ — sina ~ 1. The lightest Higgs 
boson coupling is hb — hbsin/?(l + 0{Amb/ tan^/?)), while, as happens for vanishing values 
of the bottom mass corrections, the coupling of the heaviest CP-even Higgs boson may 
become highly non-standard. 

As discussed in Sec. 2, in the large tan/? regime, the off-diagonal elements of the 
mass matrix can receive large radiative corrections with respect to the tree-level value, 
(-^ 12 ) ' ~ —{rriA + Af|)/tan/?. When both A{mb) and A{mr) are small, the coupling of 

the standard model-like Higgs boson to bottom quarks and r leptons vanishes for vanishing 
sin2Q! (see Eq. ( p.4D ). The reason for this cancellation when sin2Q; = 0 is that the standard 
model-like Higgs boson becomes a pure state, which does not couple to bottom quarks 
and r-leptons at tree level. If the bottom and r mass corrections are large, however, the 
bottom and r couplings do not cancel for sin2Q; = 0, but are just given by Ahb and Ahr, 
respectively. Indeed, from Eq. (|3.11|) (Eq. (|3.12|) ), in the limit sin a = 0 (cos a = 0), the 
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bottom coupling is given by 


r /7 \ 

hbyhij) —: ^ X 


^ Ah,. 


(3.13) 


sin/?n (1 + A(mb)) 

In this limit, the coupling to bottom quarks is much smaller than the standard model 
coupling only if |A(mb)| -Cl. A similar expression to Eq. (|3.13|) holds for the r lepton 
coupling. 

For values of A(mft) of order 1, however, a strong suppression of the bottom coupling hf, 
can still occur for slightly different values of the Higgs mixing angle a, namely 

A{mb) 


tana ~ 


tan (3 


(3.14) 


Under these conditions. 


- /A(m^) - A(mft)\ - 

'/qI IiA/^ ^ )’ 

vsmp\ l + A(mT-J J 


(3.16) 


A similar expression is obtained for the coupling h^- in the case hb = 0. Hence, if tan/3 is 
very large and A(mft) is of order one, the r Yukawa coupling may not be strongly suppressed 
with respect to the standard model case and can provide the dominant decay mode for a 
standard model-like Higgs boson. 

To recapitulate, the cancellation in the off-diagonal elements of the mass matrix can lead 
to a strong suppression of the standard model-like Higgs boson coupling to bottom quarks 
and r leptons. In general, this implies a sharp increase of the branching ratio of the decay of 
this Higgs into gauge bosons and charm quarks. However, for very large values of tan (3 and 
values of the bottom mass corrections A(mft) of order one, the branching ratio of the decay 
into T leptons may increase in the regions in which the bottom quark decays are strongly 
suppressed. 


IV. HIGGS PHENOMENOLOGY WITH LARGE tan/3 CORRECTIONS. 

A. W4) process 

The hnite corrections to the bottom Yukawa coupling are important in dehning the exact 
regions for which the bottom Yukawa coupling is suppressed. Depending on the sign of 
the bottom mass corrections and on the specihc region of supersymmetric parameter space, 
important increases or decreases in coverage may occur with respect to the case of A{mb) = 0. 
For large the coupling of the lightest CP-even Higgs boson is only slightly affected by 
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the presence of A{mb), and these corrections will not affect the discovery potential. The 
only exception is when the negative contributions to the Ai 22 matrix elements, proportional 
to hf, become relevant (see Eq. ([2.3|)). For low values of uia, instead, the bottom mass 
corrections might have an important impact in the discovery and exclusion reach for a given 
choice of parameters. 

Figures show the impact of the bottom mass corrections on the discovery reach of 
the CP-even Higgs bosons in the Wcj) channel for the case of maximal mixing, p, > 0, and 
both signs of the bottom mass corrections, assumed to be given by Ami, = K x tan/?, with 
K = ±0.05 and ±0.01. Since the Higgs sector parameters depend only on the size of the 
mixing parameters and on the sign of /i x At, while the bottom mass corrections depend also 
on the sign of /i x Mg, one can have either sign for the bottom mass corrections, for fixed 
values of the stop mixing parameters. Although the most generic features of the discovery 
reach plots are not changed by the presence of the bottom mass corrections, positive bottom 
mass corrections tend to reduce the bottom Yukawa couplings and increase the luminosity 
needed for a Higgs boson discovery. The opposite happens in the presence of negative mass 
corrections. Observe that, for values of K = —0.01, there is an improvement of the discovery 
reach at very large values of tan /?. This improvement is related to a decrease in the lightest 
CP-even Higgs boson mass induced by the negative corrections to the ^4^2 Higgs mass 
squared matrix elements, proportional to h'^, which become enhanced for large values of 
tan f3 and negative values of K. As we shall discuss below, in minimal supersymmetry 
breaking models, the bottom mass corrections tend to be positive, and hence the reach of 
the Tevatron is negatively affected. The same plots for the case of minimal mixing do not 
show as much sensitivity. 

The reason why positive bottom mass corrections suppress the reach of the Tevatron 
collider can be easily understood by studying the behavior of the effective bottom Yukawa 
couplings, Eqs. (|3.11|) and (|3.12|) . Indeed, for mA -C m^sM, since tan a/tan/? < 0, the 
expression between parenthesis in Eq. ( 3.12 ) is positive. For a hxed value of the angles 
a and /?, a positive A{mi,) tends to reduce the value of h},- For ^ m^^sM, instead, 
since tana x tan/? < 0, the effect of the bottom mass corrections on the value of hb, Eq. 
( fj.ll| ), depends on whether | tana x tan/?| is greater or less than 1. In the cases displayed 
in Figs. 9-12, this factor is always larger than one and a positive bottom mass correction 
reduces the value of hb- 


19 











B. bb(t) 


The Yukawa coupling corrections discussed above affect the associated production of a 
Higgs boson with b quarks, where the Higgs boson subsequently decays to a heavy flavor 
dual state. Indeed, the cross section times branching ratio for this process at an hih 2 hadron 
collider satisfies 


a{hih 2 bb(j){-^ bb) + X) oc hi x Bb, 


(4.1) 


and 


a{hih 2 bb(j){—> ) + X) oc hi x B^-, 


(4.2) 


where hb = hb^, hb, or hb depending on cj), and Bb and B^ are the corresponding branching 
ratios of the 0 decay into bottom quarks and r leptons, which are computed using the 
modified couplings hb and hr- In general, while the four b final state, Eq. (B. is strongly 
affected by A{mb), the bbT^T~ final state [^], Eq. ([4.2|) , is only mildly affected due to a 
cancellation of the dependence of the production cross section times branching ratios on this 
factor. 


1. On the CP-even Higgs boson masses at large values o/tan /3 


One interesting feature of the large tan (3 regime is that the CP-odd and one of the two 
CP-even Higgs bosons have similar masses and couplings. One might be tempted to take the 
signal from bbA production and decay and double it to account for the other non-SM-like 
Higgs boson. However, this approximation is optimistic, and not necessary. For example, 
when both |/i| and |a| > (P(l),this might be a poor approximation [^]. Indeed, the CP-even 
Higgs boson with similar properties to the CP-odd one has a mass approximately equal to 
Afii, with the form 


Ml, ~ 




m. 






(4.3) 


where we have omitted the two-loop corrections, Eq. (|2.3|) . A particularly interesting case 
to analyze is the maximal mixing case, a3 = 6, when the radiative corrections to ml^sM are 
maximized and Ad 12 receives only small radiative corrections for moderate or large values 
of rriA- For small mass differences compared to the average mass, one gets approximately 


rriA-mhCoi 


3mf fA 

2mA 


(a^ = 6 ). 


(4.4) 
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For \fi\ > 1, both CP-even Higgs boson masses can be significantly different from the CP- 
odd one. Fig. shows the minimal and maximal mass difference of the CP-even Higgs 
boson mass nih ijnu) with the CP-odd Higgs boson mass, for values of sin^(/3 — a) < 0.5 
(sin^(/? — a) > 0.5), in the maximal mixing case and \fi\ = 1 . A scan was performed 
over values of tua > 80 GeV. As is clear from the above expression, the maximal and 
minimal mass differences are obtained for the minimal and maximal values of chosen. 
For instance, for values of ttia — 80 GeV, d = a/O and p = ±1, one obtains a mass difference 
of about 5 GeV for large tan/3, which coincides with the results presented in the figure. 
Had we scanned over lower values of rriA, the mass difference would have increased. In our 
analysis we consider the separate signals from A and the GP-even like Higgs boson with 
similar masses and couplings as the GP-odd Higgs boson. 


2. Simulation of bb{4> —> bb) signal and backgrounds 


The numerical results for the bbcf process are based on the study of a generic neutral 
Higgs cj) (with production and decay properties of the GP-odd Higgs boson) .0 First, we 
consider the four 6 -quark final state from the decay 0(—> bb). We performed a parton-level 
simulation based on Madgraph matrix elements for the processes pp bb(j){^ bb) + X, 
—>■ bbZ{-^ bb) + X, and —bbbb (QGD). All matrix elements are evaluated at leading- 
order, using leading-order leading-order parton distribution functions (GTEQ3L), and 
a common scale ^7^/2, where is the partonic center-of-mass energy.^ The Higgs boson 
and Z boson resonances are treated in the narrow width approximation. Next-to-leading 
order QGD corrections to these processes are expected to be important, but they have not 
yet been calculated. Since the bb production cross section at the Tevatron at NLO is almost 
doubled from the LO result, as a first estimate the signal and backgrounds in this study 


are multiplied by a factor of 2 [^. This assumption will need to be considered in more 
detail elsewhere. When Gaussian statistics apply, this increases the significance i^) of a 
signal by \/2. We assume that four 6 -quark tags are necessary to reduce backgrounds (so 


®This process was first considered at the Tevatron based on a 3 jet analysis Later, it was 


reconsidered based on a 4 jet analysis [^|. The modified results of Ref. presented in Ref. |^] 
are now in general agreement with our results when comparable. 


®The scale dependence of these results is estimated to be 25% |^|. In the same study, the reach 
of a 3 or 4 6-tag analysis is estimated to be the same. 


21 









we do not consider Z{^ bb)jj or bbjj backgrounds), and that the four-tag efficiency can 
be described by an overall factor of (.45)^ ~ .2.|^ We also assume that the efficiency for 
triggering on a four b hnal state is unity. 

The signal is dehned by the following cuts: 

• A b partons with pxi > 20 GeV, i=l,4, and \ri^\ < 2, 

• max(pTi) > 40 GeV 

• Rij > 0.8, where Rij = ^(0* - + {pi - 

For events that satisfy these cuts, the distribution of all invariant mass combinations rriij = 
\J(Pi + PjY is constructed. We use a mass resolution of <Jm/M = .08 above 100 GeV, and 
aM = ^ GeV below. We use the procedure outlined earlier to combine the signals from two 
Higgs bosons which are close in mass. 

For all Higgs boson masses considered in this study, the QGD production of four 6 -quarks 
is the dominant background. The cross section for this process is proportional to a^, and, 
hence, is sensitive to the choice of scale. Therefore, an absolute prediction of the event 
rate after cuts has a large uncertainty. By studying all rriij combinations, it is possible to 
dehne a smooth background distribution and determine the overall normalization from the 
data using sidebands. If, instead, one chooses the rriij combination closest to a hypothesized 
Higgs boson mass, then a distribution similar to the signal is sculpted from the background, 
and it becomes problematic to assess the signihcance of a mass peak. 

There may be optimal cuts to increase the signihcance for heavier Higgs boson masses. 
For rrij, > 100 GeV, we include the additional requirement that | cos 0*1 < 0.8, where 9* is 
the polar angle distribution in the rest frame of the bb pair that best reconstructs to m^. 


3. Results 


The 95% G.L. exclusion and 5a discovery potentials of the bbbb channel are illustrated 


in Fig. for different total integrated luminosities and for the case of minimal mixing. 
These results imply that a GP-odd Higgs boson (and its partner with similar properties) 
can be discovered with 30 fb“^ of data at large tan/3 if rriA ~ 200 GeV. For the same 


^^This estimate is based on the double tagging efficiency of Ref. |^. The actual efficiency will 
require a detailed analysis. 
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integrated luminosity, an exclusion contour can cover from niA — 80 GeV and tan/3 ~ 15 
up to niA — 300 GeV and large tan/3 ~ 50. However, these conclusions assume vanishing 
SUSY corrections to the bottom mass (A(mft) = 0). Fig. p!5| shows the sensitivity of these 
results (for maximal mixing and /r > 0) to SUSY corrections at large tan/3. The lines in 
Fig. 1^ show the variation of the 5a discovery contour with 30 fb“^ for Arub = K x tan/3, 
with K = ±.005 and ±.01. There is a similar variation in the exclusion contours. Glearly, 
it is difficult to make a dehnitive statement about the reach of the Tevatron in the bbbb 
hnal state with limited knowledge of the size of the bottom mass corrections, A(mfe), which 
depend on the sparticle masses and mixings. It is important to realize that, for negative 
values of K and large values of tan /3, the cross section increases due to a large increase in 
the bottom Yukawa coupling hb, but this Yukawa coupling may become too large to make 
a perturbative analysis possible. This would happen, for instance for values of \K\ > 0.015 
and tan /3 50. 

One of the generic conclusions of this study is that, in the large tan /3 regime, this process 
may be useful to test regions of parameter space which will remain otherwise uncovered. In 
general, however, the reach in the m^-tan /3 plane is reduced to a relatively small region for 
values of niA of the order of the weak scale. The exact discovery and exclusion potentials 
depend strongly on the hnite SUSY corrections to the bottom mass, which can be very large. 

4- Simulation of bh4){-^ rr) signal and backgrounds 

We have also considered the possibility of detecting the Higgs boson decays into r+r“.|^ 
Since A(mT-) is expected to be small, the bbT^T~ channel is not very sensitive to SUSY- 
induced, large tan /3 corrections. This follows from Eq. (^4.2|) in the limit that the total width 
of the Higgs boson is dominated by the bb partial width. We present a preliminary study 
here, but more work needs to be done to understand the feasibility of this channel. With 
our present understanding. Fig. shows that for < 120 GeV, the reach of this channel 
is comparable to or even slightly better than the bbbb channel, and one may expect much 
room for improvement. 

To study this process at the Tevatron, we considered only the ri(—> f + X) and r 2 (— 
i + X) decays of the T 1 T 2 pair, where £ = e or /x. This combination yields a triggerable 


11 


This channel was previously considered for Run I 


However, due to a numerical error, the 


reach in tan/3 was largely overestimated in that work |^]. 
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lepton, a narrow jet, and two 6-quarks in the final state. It also has the largest branching 
ratio. The physics backgrounds are assumed to be hhZ{-^ tt) and tt production. The signal 
bb(j){^ tt) and the background bbZ{^ tt) were simulated in a similar manner as before 
and increased by a factor of 2 to account for higher-order corrections. In addition, the 
r polarization information was included for all r decays through the Tauola Monte Carlo 
program The ti background was simulated using Pythia 6.1 with the default 

settings, and forcing W decays into e, /i or r. For this analysis, the 6-parton is treated 
as a 6-jet, but the r-jet is constructed from final state particles (vr, K, etc.). The IpT is 
constructed from the real neutrinos from W boson and r decays. 

Because the number of backgrounds is smaller than the previous case, and the signal 
and Z background have similar characteristics, it is assumed that the acceptance cuts can 
be looser. The basic cuts are: 


• 2 6 partons with p\ > 10 GeV, \rf\ < 2 . 

• 1 e or p with > 10 GeV, |p^| < 2 . 

• 1 r-jet with plp > 15 GeV, |p-^ | < 2 . 

• Rij > 0 . 6 , where i,j sum over 6 ’s, £ and r-jet. 

After these cuts, the tt events produce the largest background. However, the jets and 
leptons from t-quark decays are much harder and produce much more than the typical 
signal event. The further cuts < 60 GeV and I^t < 80 GeV are imposed to reduce this 
background. For the final numbers, the GDF r-jet reconstruction efficiency ranging from 
approximately .3 to .6 is used [0, as well as a double 6 -tag efficiency of .45 and a triggering 
efficiency of unity. 

The signal is defined by a simple counting experiment, without reference to a mass 
window. Several possible improvements could greatly increase the potential of the 66 rr 
signal. First, with adequate Ifix resolution, a mass peak can be partially reconstructed. 
This could distinguish the signal from the background for S> Mz- Second, the second 
largest branching ratio for the decay of a r pair is when both r’s decay to jets. While this 
channel would greatly enhance the signal, it requires a detailed background and triggering 
analysis beyond the scope of this work. 
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V. CONSTRAINTS FROM THE DECAY B 


Si 


As shown above, the couplings of the CP-even Higgs bosons to bottom quarks depend 
strongly on In the MSSM framework, positive or negative corrections are possible. 

However, in some specihc models, the sign of the correction is correlated with the supersym¬ 
metric contribution to the amplitude of the decay process B (6 —> sy), which, at large values 
of tan/3, is proportional to x /itan/3 This is the case, for instance, in the minimal 
supergravity model, with unihcation of the three gaugino masses. For moderate and large 
values of tan [3 


At~ ^ - 1.5Mi/2, 


(5.1) 


where Aq and M 1/2 are the boundary conditions for At and the gaugino masses, respectively, 
at the Grand Unihed Theory (GUT) scale. Unless Aq > 6 M 1 / 2 , we have At x Mg < 0. More¬ 
over, it has been shown m that unless Aq and rriQ are much larger than M 1 / 2 , the expression 
for A{mb) in Eq. (|3.4|) is dominated by the gluino contributions, which are proportional to 
p X Mg. The top squark-induced corrections, proportional to the trilinear parameter At, are 
smaller than the gluino-induced ones and tend to reduce the total bottom mass corrections. 
Hence, the sign of the bottom mass corrections is determined by the gluino-sbottom loop 
contribution, which is opposite in sign to the chargino-stop corrections to the b —>■ sj decay 
rate. Gancellation of the positive contribution of the charged Higgs boson H'^ to B (6 —»• sy) 
requires AtX fi < 0, so that in these models the bottom mass corrections A{mfy) > 0. Positive 
corrections to the bottom mass (A(mb) > 0) reduce the effective bottom Yukawa coupling 
with respect to the tree level value, Eq. (E3). which reduces the discovery and exclusion 


potential of a Higgs boson in the bbbb hnal state at the Tevatron collider (see Fig. ^). As 
discussed above, positive mass corrections have also a negative effect on the exclusion and 
discovery potential of a GP-even Higgs boson in the Wcp channel (see Figs. 0 and 1 ^) . 

In general, in the absence of flavor violating couplings of the down squarks to gluinos, the 
H(6 —>■ sy) constraint on the possible values of rriA and the stop mass parameters becomes 
strong for large values of tan 13. For low values of the GP-odd Higgs boson mass, positive 
values of /i X At are disfavored by the data. Even for negative values of p. x A^, when 
niA — Mz, the suppression induced by the chargino-stop contributions tends to be too 
small to cancel the large charged Higgs boson enhancement, unless |At x p| x tan/3 becomes 
large. This cannot be achieved by pushing the value of the p parameter to values larger 
than Ms, since this would increase the chargino masses and lower the loop effect. Hence, 
large values of tan/3 and \At\ are preferred. 
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The above mentioned constraints on the stop mass parameters for low values of niA 
can be avoided in the presence of a non-trivial down squark flavor mixing. In particular, 
non-negligible mixing parameters between the second and third generation of down squarks 
can contribute to the 6 —>• sq decay rate via gluino-squark loop induced processes. For 
instance, if the gluino contributions were the only ones leading to the 6 —>■ sq decay rate, 
the branching ratio would be given by 




B{b^ sq) ~ (x) 


Ad 

^23 


Ml 


(5.2) 


where A23 is the value of the off-diagonal sbottom-sstrange left-right term in the down 
squark squared mass matrix (which we assume to be equal to the right-left one), tb — 
1.5 X X = Mg/Ms and 

1 -I- 4a; — 5a;^ -|- 4a; ln(a;) -|- 2x^ ln(a;) \ 


F{x) = 4 


(5.3) 


Y 8(1 — xY ) 

This expression ignores the potentially relevant contributions coming from a left-left down 


squark mixing Observe that, for Mg = 0{Ms) and Ms = 0(1 TeV), the contribution 

of the gluino mediated diagram to the b ^ s'j branching ratio is of the order of {A^^/M^Y. 
Hence, even a small left-right mixing, A23, of order of 10“^ x M| can induce important 
corrections to the amplitude of this decay rate. It is straightforward to show that these low 
values of A23 do not have an immediate impact on the Higgs boson sector. 

In the presence of non-trivial flavor mixing in the down squark sector, large corrections 
to the amplitude of the 6 —^ sq decay rate may be induced. These corrections may be 
helpful in determining values of B{b ^ s'j) consistent with experimental data for small 
values of rriA and/or positive values of p x At. For the above reasons, in our presentation, 
we have decided to keep the results for both signs of At x Jl. The reader must keep in 
mind, however, that positive values of this parameter for low values of ruA would imply 
a more complicated flavor structure that the one appearing in minimal gauge mediation 
or supergravity schemes. Observe that the contributions of the charged Higgs boson and 
top-squark loops to 6 —>• sq, should be computed including the effect of the bottom mass 
corrections A{mA in the dehnition of the bottom Yukawa coupling hi,. A next-to-leading- 
order SUSY QCD calculation of B(6 —sq) can be found in Ref. . 


VI. CONCLUSIONS 


We have presented a study of some of the MSSM Higgs boson signatures of relevance 
to the Tevatron collider. We hrst analyzed the possibility of hnding the lightest or heaviest 
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CP-even Higgs bosons in the W(j) channel. Quite generally, for moderate and large tan/3, 
either the lightest or the heaviest CP-even Higgs boson has SM-like couplings to the vector 
bosons. Therefore, most of the tan (3 — rriA plane is covered in this channel by producing the 
corresponding Higgs boson with mass below 130 GeV, provided there is sufficient integrated 
luminosity, of the order of 30 fb“^. However for ruA = m^sM ~ 110 — 130 GeV and large 
tan f3, neither h nor H has SM-like couplings to the vector gauge bosons and the coverage is 
decreased. This problem is more pronounced for large values of the stop masses and mixing 
parameters. For smaller values of tan/3 ~ 5, this problematic region extends to smaller 
values of m^. In these cases, very large luminosity, above 30 fb“\ will be needed or the 
contribution of other production processes will be necessary to assure full coverage. To cover 
the region of tan/3 > 10 and ttia — m^sM, we have combined the two CP-even Higgs boson 
signals when their masses are close to each other. Another possibility may be to explore the 
hA and HA channels, which will suffer from similar suppression factors in the production 
cross sections (sin^(/3 — a) ~ cos^(/3 — a) ~ 0(0.5)), but may be combined with the Wcf) 
production process. 

Furthermore, the branching ratio for the decay of the CP-even Higgs bosons into bottom 
quarks can be very different from the standard model one. In particular, this takes place for 
moderate and large tan /3, when the off-diagonal elements of the Higgs boson mass matrix 
can be strongly modihed by radiative corrections induced mainly by top squark-loops. We 
derived an approximate expression in terms of the MSSM parameters to clarify when this 
occurs, and provided examples when the decays into bottom quarks are suppressed in a large 
region of parameter space, thereby negatively affecting the Tevatron reach. 

We have also emphasized that, due to supersymmetry breaking effects, the values of the 
bottom and r Yukawa couplings to the CP-even Higgs bosons may be different from the ones 
computed including only standard QCD corrections. Indeed, non-decoupling effects induced 
by supersymmetry breaking can become particularly important for large tan/3, leading to 
modifications in the Higgs boson discovery and exclusion potentials at the Tevatron. For 
instance, if the bottom mass corrections A(mft) are of order one, the decay of the standard 
model-like Higgs boson into may be enhanced, while decays to bb are suppressed. 

These supersymmetry breaking effects on the bottom and the r Yukawa couplings can also 
have an impact in the phenomenology of the charged Higgs boson. In particular, they can 
be relevant for determining the Tevatron limits on top decays into charged Higgs bosons at 
large tan /3 . 

The bottom Yukawa coupling corrections are particularly important for the bb(j) process, 
because the production cross section is proportional to the square of the bottom Yukawa 
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coupling. We performed a phenomenological study to investigate the relevance of these 
corrections. Even with large luminosity factors, of the order of 30 fb“^, and negative bottom 
mass corrections, which enhance the production rate, the Tevatron can discover a CP-odd 
Higgs boson (together with a CP-even Higgs boson with mass and couplings similar to 
it) only if its mass is not larger than about 200-300 GeV. The reach is only efficient for 
moderate or large values of tan (3. Such values for the CP-odd Higgs boson mass give positive 
contributions to B{b —>■ sy), and the discovery of such a Higgs boson would constrain the 
masses and mixing angles of the top squarks unless a non-trivial mixing between the second 
and third generation down squarks is present. 

The computation of the Higgs boson mass matrix elements considered in this article |§] 
is still affected by theoretical uncertainties, most notably, those associated with the two- 
loop, hnite, threshold corrections to the effective quartic couplings of the Higgs potential. 
Recently, a partial, diagrammatic, two-loop computation of the Higgs mass has been per¬ 
formed 1^. In the limit of large tua, these additional contributions lead to a slight mod- 
ihcation of the dependence of the lightest CP-even Higgs boson mass on the stop mixing 
parameters. For instance, although the upper bound on the lightest CP-even Higgs boson 
mass for squark masses of approximately 1 TeV is approximately the same as the one ob¬ 
tained to next-to-leading-order accuracy (as done in this work), the upper bound on the 
Higgs boson mass is reached for values of |a| ~ 2 instead of |a| = \/6, and has a weak 
dependence on the sign of d. A diagrammatic computation of the two-loop corrections in¬ 
duced by the top Yukawa coupling, which are included at the leading-logarithmic level in 
our computation, is, however, still lacking. 


Summarizing, at present, the Wcj) channel with the W decaying leptonically and the 
Higgs boson decaying into b quarks remains the golden mode to test the MSSM Higgs sector 
at the Tevatron. The other channel we have analyzed, bb(j) production with the subsequent 
decay of 0 into b quarks and r leptons, proves to be very useful to cover regions of large 
tan (3 and small to moderate up to about 250 GeV. However, the reach in these channels 
requires a large total integrated luminosity. Because of this, other production processes and 
Higgs boson decay modes need to be carefully investigated if we want to fully probe the 
MSSM Higgs sector at the Tevatron. We have identihed regions of parameter space where 
the Higgs decay into bb is strongly suppressed. In these regions, other search techniques 
will have to be considered, due to the presence of enhanced decays to W*W*, gg and cc 
hnal states. Glearly, there is a motivation to study these hnal states, and in particular 
the W*W* one, even for lighter Higgs bosons for which the SM Higgs boson decay rate is 
strongly suppressed. In addition other production processes like the associated production 





of hA oi HA with the subsequent decays into b quarks and r leptons may also be useful. 

A careful study of all different possibilities, which may be relevant in different regions of 
parameter space, and the combination of channels may allow a full coverage of the MSSM 
parameter space with luminosities achievable at the Tevatron. If that is the case, the Teva- 
tron can discover a light Higgs boson which might be beyond the presently expected LEP2 
reach for generic values of the supersymmetric mass parameters. Most importantly, the 
detection of one or more Higgs bosons at the Tevatron will give very valuable information 
about the Higgs and stop sectors of the MSSM. 

In the hnal stages of this work, two preprints appeared on related topics. One addressed 
the issue of the Wcp reach of the Tevatron collider The other commented on the possible 
effects of large tan/3 corrections to bbcj) production at hadron colliders [31|. In the special 
cases when the analyses are comparable we tend to agree with their results, although the 
authors of Ref. do not see any visible dependence on the sign of /i x a. The present 
work goes beyond those studies by providing a detailed numerical and theoretical analysis 
of the dependence of the Tevatron discovery potential on the MSSM parameter space. 
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FIGURES 


V(|) Production 

V=W at the Tevatron, V=Z at LEP2 



FIG. 1. 95% C.L. bound on i? = ^ function of the Higgs boson mass 

for the Tevatron and LEP2. 
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c(W^)BR(^^bB)/c(W^3M)BR(^SM^bB) 


V(|) Production 

V=W at the Tevatron, V=Z at LEP2 



(GeV) 


FIG. 2. Same as Fig. 1 except for 5 cj discovery. 






CQ. 
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Channel, M 3 =l TeV, |a<0 

8 fb'^ 16fb'^ ^ 24fb‘^ 32 fb'^ 



100 150 200 250 300 350 400 

(GeV) 

FIG. 3. 5(T discovery contours for the IF0(—> bb) mode at the Tevatron in the MSSM for 
maximal mixing, // < 0, and Ms=l TeV. Different shadings correspond to different integrated 
luminosities. 
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Channel, M 3 =l TeV, |a>0 
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FIG. 4. Same as Fig. 3 but with fi > 0. 
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FIG. 5. Same as Fig. 3 but for minimal mixing. 
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Channel, A^.^, M^=l TeV 
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FIG. 6. Same as Fig. 5 but using the present CDF mass resolution. 
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Channel, A=-|a=l TeV 


8 fb'^ 16fb'^ ^ 24 fb"^ 32 fb'^ 
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FIG. 7. 5a discovery contours for the VF(^(—> bb) mode at the Tevatron for A = —/r = 1 TeV 
and Ms = 1 TeV. 
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Channel, A=-|a=1.5 TeV 


8fb'^ 16fb'^ ^ 24 fb"^ 32 fb'^ 



100 150 200 250 300 350 400 

(GeV) 

FIG. 8 . Same as Fig. ^ but for A = —/i = 1.5 TeV. 
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W(|) Channel, Am^=-.005 tan^, |a>0 


8 fb'^ 16fb'^ ^ 24 fb"^ 32 fb'^ 
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(GeV) _ 

FIG. 9. 5(7 discovery contours for the W(f){-^ bb) mode at the Tevatron for maximal mixing and 
fj. > 0 after including large tan/3 corrections: A{mb) = —.005tan/3. Different shadings correspond 
to different integrated luminosities. 
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W(|) Channel, Am^=-.01 tan^, |a>0 
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FIG. 10. Same as Fig. ^ but for A(mft) = —.Oltan/3. 
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W(|) Channel, Am^=.005 tan^, |i>0 
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FIG. 11. Same as Fig. ^ but for A{mb) = .005 tan/I. 
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W(|) Channel, Am^=.01 tan^, |a>0 
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FIG. 12. Same as Fig. ^ but for A(mft) = .01 tan/3. 
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Bounds on the Higgs Boson Mass Difference 



FIG. 13. Mass difference between the CP-odd Higgs boson and the CP-even Higgs boson (with 
properties similar to the CP-odd one) as a function of tan f3 for several MSSM parameter choices. 
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FIG. 14. 95% C.L. exclusion and 5a discovery contours for bbcj) production at the Tevatron in 
the MSSM for minimal mixing. The shaded areas correspond to different exclusion and discovery 
contours at different integrated luminosities for the bbbb final state. The different lines show the 
same for the bbT~^T~ final state, except that the lower solid line shows the 95% C.L. exclusion 
contour for the bbbb final state and 30 fb“^. 


Limits on bD(|) Production 
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Large tan|3 Corrections to bBbB Limits 



(GeV) 

FIG. 15. 5a discovery contours for the bb(j){-^ bb) mode and 30 fb“^ with maximal mixing, 
/i > 0, and A(mfe) = (0, ±.005, ±.01) tan/3. 
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